arXiv:1506.08748v2 [physics.ins-det] 15 Sep 2015 


Preprint typeset in JINST style - HYPER VERSION 


Measurements of wavelength-dependent double 
photoelectron emission from single photons in 
VUV-sensitive photomultiplier tubes 


C.H. Faham", V.M. Gehman"* A. Currie"^, A. DobP, P. Sorensen" and R.J. GaitskelP 

“Lawrence Berkeley National Laboratory, 1 Cyclotron Rd., Berkeley CA 94720, USA 
^Imperial College London, High Energy Physics, Blackett Laboratory, London SW7 2BZ, UK 
“Brown University, Dept, of Physics, 182 Hope St., Providence RI02912, USA 


Abstract: Measurements of double photoelectron emission (DPE) probabilities as a function 
of wavelength are reported for Hamamatsu R8778, R8520, and R11410 VUV-sensitive 
photomultiplier tubes (PMTs). In DPE, a single photon strikes the PMT photocathode and produces 
two photoelectrons instead of a single one. It was found that the fraction of detected photons that 
result in DPE emission is a function of the incident photon wavelength, and manifests itself below 
~250 nm. Eor the xenon scintillation wavelength of 175 nm, a DPE probability of 18-24% was 
measured depending on the tube and measurement method. This wavelength-dependent single 
photon response has implications for the energy calibration and photon counting of current and 
future liquid xenon detectors such as EUX, EZ, XENONIOO/IT, Panda-X and XMASS. 
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1. Introduction 

Photomultiplier tubes (PMTs) are currently one of the most widely used photodetectors due to their 
large collection area and high gain. In radiation detectors using xenon as an active medium, vacuum 
ultraviolet (VUV)-sensitive PMTs are used to detect the 175 nm xenon scintillation light OJ, 0]. For 
photon-counting and calorimetry applications, the single-photon response of PMTs must be well 
characterized to understand the overall detector behavior. 

The most important contribution to the PMT photon detection efficiency is the quantum 
efficiency (QE), which is defined as the ratio of electrons produced in the photocathode (referred to 
as photoelectrons, or phe) divided by the number of incident photons [0]. QE is often incorrectly 
treated as the probability of an incident photon to produce a photoelectron, but if more than one 
photoelectron is produced by a single photon, this equivalence is no longer valid. QE is a function 
of the incoming photon wavelength due to a number of competing factors, including the PMT 
window transmission probability and the photon absorption length in the photocathode. QE is 
typically limited to ~30% at 400 nm, with decreasing sensitivity at lower wavelengths. Eor PMTs 
specially designed to be sensitive to the VUV band, the QE shows a mai'ked enhancement below 
200 nm and above the transmission cutoff of the synthetic silica window at 160-170 nmfll], 0, 0]. 
The present work was prompted by observed differences in the PMT pulse height and area spectra 
for single photons of xenon scintillation light and blue EED light in EUX[0, 0]. Specifically, a 
feafure like the one discussed in Section 0 of this document, led us to start investigating not only 
the efficiency, but also the shape of the single photon response of EUX PMTs as a function of 
incoming wavelength. 

The process of emitting two photoelectrons from the photocathode from a single incoming 
photon, referred to as double photoelectron emission (DPE), has been studied before in the context 
of the photoelectric effect on metallic surfaces [0, EL 0]. PMT photocathodes are expected to 
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exhibit a similar behavior, provided that the incoming photon is sufficiently energetic (greater than 
twice the photocathode work function). We have found in this study that for wavelengths lower 
than ~250 nm, the DPE probability is non-negligible for all three VUV-sensitive Hamamatsu 
R11410 (3-inch circular photocathode), Hamamatsu R8778 (2-inch circular) and R8520 (1-inch 
square). The R8520 PMT was used in the XENONIO, XENONIOO, Panda-X (top array) and EZ 
(xenon veto region) experiments. The R8778 was utilized in the EUX experiment, while a variant 
(R10789 hexagonal head) was used in the XMASS experiment. Eastly, the larger R11410 has been 
used in the Panda-X experiment (bottom array) and will be used by next-generation dark matter 
experiments, such as EZ, XENONIT/NT, and RED, among others. 

2. Experimental Setup 



Figure 1. Experimental setup for the DPE measurement. 


The experimental setup used for all measurements is shown schematically in Eigure |T] and 
is similar to the one used in References HIO. ITTl. IT2I1 . A Model 632 Deuterium Eight source, from 
McPherson, Inc. 11131] created a broad spectrum extending from 110-950 nm with the window facing 
into a vacuum space pumped down to below lO^'* Torr. A focusing mirror imaged the lamp window 
on a holographic diffraction grating inside a Model 234/302VM 0.2-Meter Monochromator, also 
from McPherson, Inc. The grating separated the incoming light in reflected angle by wavelength 
and could be rotated to select a specific wavelength through the exit of the monochromator. The 
monochromator body has adjustable slits (5 mm tall by 0-3 mm wide) at its entrance and exit. To 
first order, the entrance slit only adjusts the intensity of the light exiting the monochromator, while 
the exit slit also adjusts the width in wavelength of the monochromator output. A narrower exit 
slit samples a smaller range of reflected angles from the grating and therefore results in a narrower 
wavelength distribution. Eor the duration of this study, the exit slit was set to 0.5 mm and the 
entrance ranged from 0.2-0.5 mm to regularize the single photon rate as measured by counting 
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peaks above threshold in the data stream. In this configuration the, full-width, half-maximum 
of the light source output was measured to be 4.4 nm between 200 and 300 nm with an Ocean 
Optics QE65000 UV/Vis spectrometer (which is only sensitive above 200 nm). Second-order (half¬ 
wavelength) peaks are present in the monochromator setup. Because the synthetic silica windows 
in all three PMT models are perfectly opaque to light below 155 nm, the maximum wavelength 
used in this study was 300 nm, which ensured no second-order peaks were detected by any of the 
PMTs. All of the PMT measurements were performed at room temperature. 

Table 1. Bias and acquisition settings for each of the 3 tested PMTs. 

Mean single phe 


PMT type 

height [mV] 

DAQ threshold [mV] 

1-inch R8520 

9.2 

0.75 (0.082 single phe) 

2-inch R8778 

36.8 

1.55 (0.042 single phe) 

3-inch R11410 

21.2 

1.40 (0.066 single phe) 


The PMTs examined in this study were mounted in black Delrin® holders to provide electrical 
insulation from their surrounding grounded flange, as well as to minimize reflected light. PMT 
high voltage bias and signal were fed in and read out through coaxial feedthroughs on the vacuum. 
PMT signals were digitized by a PicoScope 5000 series oscilloscope. The PicoScope was triggered 
every 10 ms by an external pulse generator, and read out using a custom MATLAB [1141] data 
acquisition (DAQ) code. When the PicoScope received a trigger signal, the DAQ code read in 
a 5-ms long buffer, which was sampled at 500 MS/s (2 ns samples), and pulses above threshold 
were selected. The pulse threshold was set depending on noise levels and gain for each PMT, and 
ranged from 0.75-1.55 mV. Threshold and single PE peak values are reported in Table |T]. Eor each 
pulse above threshold, a waveform of 150 samples (300 ns) centered around the pulse maximum 
was selected and was then baseline-corrected by subtracting the average of the first 10 samples. 
This pulse selection was performed in order of descending pulse size and the selection finished 
when there were no pulses left in the 5 ms buffer with a height above the set threshold value. 
Each acquisition (one per wavelength) was saved to a .mat file with the pulse waveform data and 
acquisition metadata. A total of 1 million pulses were collected at each wavelength for each PMT. 
The rate at which pulses were collected is plotted in Eigure @ and was calculated by dividing the 
number of observed pulses by the total amount of time searched. At these rates, the pileup fraction, 
i.e. the fraction of 300 ns windows containing >1 pulse, ranged from 0.9% to 1.5%. 

3. Results 

The waveforms were loaded and analyzed in Python using SciPy flT^II and pyROOT llTHl . Single 
photon response spectra for each PMT at each wavelength were created by histogramming the 
waveform area. We show the qualitative variation of these pulse area histograms for the R11410 
as a function of wavelength in Eigure ^ The spectra in Eigure 0 were normalized by the integral 
in the region between 50-120 mVns to better illustrate the increase of the amplitude of the double 
photoelectron peak with decreasing photon wavelength. 
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Figure 2. Rate of single photon and dark pulses above DAQ threshold during data acquisition for all 
wavelengths and PMT models (top panel). Bottom panel shows the same photon rate data as the top, but 
plotted on a linear scale for easier comparison. 


We further illustrate this difference in Figure ^ where we have plotted an overlay of the first 
2000 pulses collected at both 300 nm (left) and 175 nm (right) for the R11410. The 175 nm 
waveforms clearly have a higher amplitude and effectively the same width compared to the 300 nm 
waveforms. The photon rate recorded by our DAQ was nearly identical: 46.1 kHz at 300 nm and 
45.4 kHz at 175 nm. The central region of 0is centered at just less than 100 ns. Note that there 
are a few dozen pulses that fall outside of the central time. This is consistent with our estimate of 
the pileup rate, and these pulses were removed from our subsequent analysis with a simple pulse 
quality cut to remove events with more than one significant local maximum. 

Another feature visible in both Figures ^ and ^ is a population of pulses aiising from incident 
photons which are not absorbed by the photocathode, but create a photoelectron on the first dynode 
instead. These events populate the peak in the pulse area histogram centered around 20 mV ns 
in Figure ^ and the persistent trace visible at 1.5 mV in amplitude in Figure ^ We verified this 
by running with the photocathode and first dynode pins on the PMT base shorted together and 
the overall bias adjusted to maintain the same field befween fhe resf of the dynodes as in regular 
operating conditions. In this mode, the efficiency for collecting electrons from the photocathode is 
very small, while the efficiency for collecfing elecfrons from the first dynode remains unchanged. 
When running in this shorted photocathode configuration, the mean and width of this 20 mV ns 
distribution remain unchanged while the single and double phe distribution vanish. The first dynode 
photoelectrons should in principle, exhibit a similar DPE behavior as those from the photocathode. 
This is an interesting topic for future study, but lies outside the scope of this article. 


- 4 - 


























C/) 

' 

c 

3 

o 

o 

"D 

CD 

"S 

E 



0 100 200 300 400 500 600 

Pulse area [mVns] 


Figure 3. Superposition of pulse area spectra for the R11410 PMT for different wavelengths. Each spectrum 
is normalized by the integral in the region between 50-120 mV ns in order to show the effect more clearly. 
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Figure 4. First 2000 pulses at 300 nm (left) and 175 nm (right) for the R11410 PMT acquisition. The pulse 
rates were 46.1 kHz and 45.4 kHz, respectively. All acquired waveforms are shown, but pulse quality cuts 
to remove pileup pulses were applied for the analysis. Direct photon hits on the PMT first dynode can be 
seen at approximately 1.5 mV in height. The rates of both first dynode photoelectron emission and double 
photoelectrons, relative to that of single photoelectrons, are higher at lower wavelengths. 
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The pulse area histograms were used to calculate the fraction of PMT pulses exhibiting double 
photoelectron emission as a function of wavelength. Two different methods were employed: 

1. Perform a simultaneous single and double photoelectron distribution fit for each spectrum, 
and use the fitted areas to estimate the DPE fraction. 

2. Calculate the shift in the mean pulse area, which should directly correspond to the DPE 
fraction. 


Eor the first method, our fit model was a pair of Gaussian functions^ plus a constant 
background offset: 

~^sphe ^ ^ ^~^dphe 

“ 1 “ ^dphe^ 


F{x') Nsphe^ 


°sphe 




^^dphe j 


+ C, 


(3.1) 


where Nsphe and Ndphe are the normalizations, Psphe and Pdphe are the mean values, Osphe and Odphe 
are the widths for the single and double photoelectron peaks, and C is the constant background 
term. The initial guesses for each fit come from a Gaussian fit to the region very near the single 
phe peak in the 300 nm spectrum, which should have the smallest DPE fraction for the acquired 
data. The double phe peak parameter guesses are then calculated under the assumption that the 
double phe population is defined by sampling fhe single phe populafion fwice, i.e. Pdphe = 2 Psphe 
and Odphe = Osphe- Bofh means and widfhs were fhen allowed fo floal wifhin ±5% of fhe 
initial guess value. The DPE fraction for each specfrum is fhe area of fhe double phe disfribufion 
divided by fhe sum of fhe area of fhe single plus fhe double phofoelecfron disfribufions. We show 
an example of a til fo fhe R11410 pulse area specfrum af 175 nm in Eigure for which fhe DPE 
fraction was calculated to be 22.5 ± 1.0%. 

To calculate the DPE fraction using the shift in the mean of the pulse area spectrum, we 
truncate the low side of the spectrum at the valley between the single phe peak and the first dynode 
and pedestal structure, e.g. around 50 mVns for the spectrum in Eigure H. We then simply take the 
bin-weighted mean of each pulse area spectrum. The mean shift is calculated relative to the 300 nm 
spectrum by the formula: 


Mean Shift(A) 


Pspec (2-) Pspec (2- — 300 nm) 
Pspec (2- — 300 nm) 


(3.2) 


where Pspeci^) is the mean of the pulse area spectrum at wavelength A. This quantity is therefore 
identically zero by construction at 300 nm, and allows for a simple, model-independent estimate of 
the DPE fraction versus wavelength. 

We plot the DPE fraction using both the fit (in the left panel) and mean shift (in the right panel) 
methods for each PMT in Eigure ^ 

Uncertainties in the traces of Eigure ^ are a quadrature sum of two effects. One is the intrinsic 
uncertainty from the process of calculating the DPE fraction. This incorporates the fit uncertainties 
in the left panel and the counting statistics uncertainties for the right. The second one is the 

'while PMT pulse area distributions can have significant deviations from Gaussian in the limit of sufficiently large 
statistics, the agreement between this model and this data is excellent. The Gaussian-based model used in this article is 
much easier to interpret than others, such as Polya or scaled Poisson distributions because it fits the mean and width of 
each peak directly. 
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Figure 5. Fit to the single, and double photoelectron spectrum (top panel) and fit residuals (bottom panel) 
for the R11410 at 175 nm. Residuals are normalized to the bin-by-bin uncertainty (assuming simple Poisson 
counting statistics) and are plotted over the fit range. The DPE fraction for this fit to the spectrum is 22.5 ± 
1 . 0 %. 



Figure 6. DPE probability for all three PMTs tested as a function of incident photon wavelength using the 
fit (left) and mean shift (right) methods. 


systematic uncertainty from the range over which each method was carried out, and this is the 
dominant contribution to the overall uncertainty. For the fit, this is just the uncertainty arising from 
changes to the fit range. For the mean shift, this is the uncertainty arising from adjusting the range 
over which the histogram mean was calculated. Both methods agree in the overall DPE fracion 
vs. wavelength trend, and point to a DPE fraction that turns on between 250 and 260 nm, climbs 
steadily with decreasing wavelength, and reaches its maximum value of approximately 20% near 
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185 nm. 


4. Implications for QE Measurement and Photon Counting in Xenon Experiments 

PMTs are commonly used to detect photons in the visible band, and visible light sources are more 
readily available for bench-top measurements. This means nearly all PMT characteristics are 
measured at visible wavelengths, usually with blue (400^50 nm) light emitting diodes (LEDs). 
The fact that a single photon can cause more than one photoelectron to be emitted from a 
photocathode at VUV wavelengths makes it important to measure the relevant PMT characteristics 
in the VUV. It is therefore important to distinguish two different figures of merit in place of a single 
quantum efficiency rj. Lef be fhe ratio of phofoelectrons emiffed from fhe phofocafhode 

fo phofons incidenf, and lef t]p(A) be fhe fracfion of incidenf photons generafing one or more 
phofolecfrons. The former is fhe mean response and fhe latter is roughly (ignoring firsf dynode hifs 
for fhe momenf) fhe probabilify of a defecfable response; in fhe regime where fhe emission of fhree 
or more phofoelecfrons is negligible, fhey are relafed by: 

T]^(A) = (1 + DPE fraction) T]p(A). (4.1) 

The disfincfion befween fhese fwo efficiencies fogefher wifh fheir wavelengfh dependence have 
several implicafions for experimenfs which measure VUV scintillation lighf wifh an array of single- 
phe sensitive PMTs. 

The mean response t]p(A) can be measured via fhe phofocafhode currenf under illumination 
by a calibrated light intensity, and is usually quoted in PMT spec sheets. However, it is T]p(A) that 
often determines a detector’s energy threshold via the probability for a given primary scintillation 
signal to meet an n-PMT coincidence requirement. Simulations which use measurements of 
as if they were probabilities for detection will over-estimate efficiencies of fhis fype whenever fhe 
DPE fracfion is greater fhan zero. 

Away from fhreshold, DPE affecfs bofh resolution and fhe calibrafion of absolufe yields of 
quanfa. Eor a given mean number of phofoelecfrons, fhe resolution of short-wavelength signals is 
worsened, with respect to what would be expected from longer-wavelength calibration, by both the 
smaller underlying binomial detection probability and the higher fractional standard deviation of 
the single-photon response. We suggest that t]p(A), calculated for instance via equation p~T| , be 
quoted in the specification of VUV-sensitive photomultipliers to enable more accurate modeling 
of thresholds and resolution. Calibration of the single photon response using VUV signal light, as 
in lEIl, instead of longer wavelengths allows pulse areas to be interpreted directly as a number of 
detected signal photons without separate measurement of the DPE fraction. Simple energy scale 
calibrations that derive from the comparison of signals summed over all PMTs to mono-energetic 
calibration lines from known radioactive sources (such lines are typically much higher energy than 
the threshold) should be largely insensitive to the DPE fraction because it is in both the numerator 
and denominator of that energy calibration ratio. Energy calibrations which depend on a model of 
detector resolution, especially those using continuous spectra, must however account for the DPE 
fraction. 

The relatively large variance of the VUV photon pulse area distribution, due to DPE, argues for 
counting integer photons in waveforms, in addition to integrating pulse areas, to optimize resolution 
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of small scintillation signals. A further advantage is that first-dynode photoelectrons could be 
used to boost efficiency in the same few-photon regime. The first-dynode peak seen in Figure 0 
contributes approximately 20% of the photons above threshold at 175 nm. Further study is needed 
of how much additional light collection can thus be gained in full detectors: the present results are 
for normally incident light in vacuum and the probability of transmission through the window and 
subsequent absorption on the first dynode depends on the surrounding material index of refraction 
and incident photon trajectory; additionally, dynode voltages could be tuned to raise more of the 
first dynode signal above the photon-counting threshold. 

5. Summary and Future Work 

This study addresses a number of interesting questions that impact numerous physics programs 
that use xenon as a radiation detection medium. The normalization of PMT signals to count 
photoelectrons is core to energy calibration and acceptance characterization, and this article 
addresses a number of important wavelength dependent issues. The magnitude of the effect shown 
in Figure ^ is similar to that seen in LUX single-photon spectra and may account for the short- 
wavelength enhancement of in VUV-sensitive PMTs flU]. 

There are a number of follow-up measurements that could be done to further investigate the 
emission of two photoelectrons into the VUV. First, there is a 10-15% relative variation in the DPE 
fraction at the shortest wavelengths to which these PMTs are sensitive. Because only one PMT 
of each model was tested for this study, it remains unclear whether this corresponds to variation 
within or between PMT models. A study involving a larger number of each device would help to 
address this question. Furthermore, most xenon-filled defectors operafe predominafely in fhe liquid 
phase, meaning fhaf fhe PMTs musf operafe af cryogenic femperafures (165 K). The fesfs for fhis 
sfudy were performed af room femperafure, so a fulure sfudy examining femperafure dependence 
of fhe DPE fraclion would be very informalive. Anofher subjecf for fulure sfudy would be to check 
if fhe calibrafed phofodiodes used for absolute QE measuremenfs exhibif fhis same DPE effecl, and 
how fhaf affecfs fhe wavelengfh-dependenf QE measuremenfs lhal the manufacturer provides. 
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